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ABSTRACT 

NGC 1846 and NGC 1783 are two massive star clusters in the Large Magellanic 
Cloud, hosting both an extended main sequence turn-off and a dual clump of red 
giants. They present similar masses but differ mainly in angular size. Starting from 
their high-quality ACS data in the F435W, F555W and F814W filters, and updated 
sets of stellar evolutionary tracks, we derive their star formation rates as a function 
of age, SFR(t), by means of the classical method of CMD reconstruction which is 
usually applied to nearby galaxies. The method confirms the extended periods of 
star formation derived from previous analysis of the same data. When the analysis is 
performed for a finer resolution in age, we find clear evidence for a ~ 50-Myr long 
hiatus between the oldest peak in the SFR(i) , and a second prolonged period of star 
formation, in both clusters. For the more compact cluster NGC 1846, there seems to 
be no significant difference between the SFR(i) in the cluster centre and in an annulus 
with radii between 20 and 60" (from 4.8 to 15.4 pc). The same does not occur in 
the more extended NGC 1783 cluster, where the outer ring (between 33 and 107", 
from 8.0 to 25.9 pc) is found to be slightly younger than the centre. We also explore 
the best-fitting slope of the present-day mass function and binary fraction for the 
different cluster regions, finding hints of a varying mass function between centre and 
outer ring in NGC 1783. These findings are discussed within the present scenarios for 
the formation of clusters with multiple turn-offs. 

Key words: Stars: evolution - Hertzsprung-Russcll (HR) and C-M diagrams 



1 INTRODUCTION 

The star clusters NGC 1846 and NGC 1783 represent pro- 
totypes of massive intermediate-age star clusters in the 
LMC containing multiple main sequence turn-offs (MM- 
STO; IMackev fc Brobv Nielsen! 120071; IMackev et all |200S| ; 
iMilone et all |2009t IGoudfrooii et al.l l2009h . They have 
masses of about 1.5 and 1.7 x 10 s M@, respectively, which 
locates them among the most massive LMC clusters excep t 
for the old globulars (see e.g. fig. 13 in lGirardi et al.ll 19951 ). 
In addition to the MMSTOs, they also seem to present a 



dual clum p of red giants, in similarity to the SMC cluster 
NGC 419 (|Girardi. Rubele fc Kerberl(2009h and the LMC's 
NGC 1751 (|Rubele et al.ll201ll ). 

The presence of MMSTOs is commonly inter- 
preted as the signature of continued star formation, 
or multiple events of st ar formation, spanning a few 
100 Myr in ti me (e .g. IMackev fc Brobv Nielsen! |2007| ; 
Mackev et ail 120081; IGoudfrooii et al.l |2009[ 12011a 



Conrov fc Spergell 1201 it Girard i. Eggcn berger fc Miglid 
201 it iKeller. Mackev fc Da Costal l201ll fl. For a well- 



* Based on observations with the NASA/ESA Hubble Space Tele- 
scope, obtained at the Space Telescope Science Institute, which 
is operated by the Association of Universities for Research in As- 
tronomy, Inc., under NASA contract NAS5-26555 



1 The disp ersion in rotational veloci ties in a coeval cluster, ad- 
vocated bv lBastian fc de Minlj d2009h , was shown by not to pro- 
duce MMSTOs similar to the observ ed ones, both theoretically 
dGirardi. Eggenberger fc Miglid l201ll ) and observationally (see 
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defined age range between 1.2 and 1.7 Gyr, star clusters 
with a turn-off age spanning a few 100 Myr will temporarily 
contain stars that ignited helium under both non-degenerate 
and degenerate conditions, and hence naturally develop a 
dual red clump (jGirardi. Rubele fc Ker"be3l2009h . 

The main difficulties with the prolonged-star forma- 
tion history (SFH) interpretation of MMSTOs, are related 
with the theories of star formation and gas dynamics in- 
side the relatively shallo w po tential wells of s t ar clus ters. 
IConrov fc Spergell (|201lh and iGoudfrooii et al.l (|2011ah de- 
scribe scenarios for the continuation of star formation over 
long timescales which impl y that more massiv e clust ers may 
have more extended SFHs. IConrov fc Spergell (|201ll ) specif- 
ically note that all clusters wit h MMSTOs have masses 
higher than ~ 10 4 Mq, whereas IGoudfrooii et al.l (|2011al l 
find a correlation between the estimated escape veloci- 
ties at an age of 10 Myr and the concentration of stars 
in the brightest half of th e MM STO region. Moreover, 
iKeller. Mackev fc Da Costal (|201ll ) notice the correlation 
between the cluster core radius r c and their SFHs, in the 
sense that all known clusters with MMSTOs have r c > 
3.7 pc. Further testing these trends, and finding additional 
correlations with other cluster parameters, are probably nec- 
essary steps to clarify the origin of clusters with MMSTOs. 

In this context, the pair of clusters NGC 1846 and 
NGC 1783 is extremely interesting. While they have 
very similar total masses ( log(M/M a ) = 5.17 ± 0.09 
and 5.25 ± 0.09, respectively; IGoudfrooii et al1l2011bh and 
mean ages (1.73 ± 0.10 and 1.70 ± 0.10 Gyr, respectively; 
IGoudfrooii et al.ll2011aT ). and a similar location in the North- 
west portion of the LMC (hence similar fore/background, 
and likely the same distance), they have very different an- 
gular sizes: for NGC 1846 the core radius is r c = 26.0" = 
6.3 pc, and the conc entration index c = r t /r c (where r t is the 
tidal radius) is 6.2 jGoudfrooii et al.ll2009l ); for NGC 1783 
these quantities are r c = 37.7" = 9.1 pc and c = 9.2 
IGoudfrooii et ai1l2011bl ). respectively. So, in these clus- 
ters we can test whether there is any measurable difference 
in their SFHs that can be interpreted as a result of the 
different radii, and in the light of the correlation noted by 
IKeller. Mackev fc Da Costal (|201ll ). 

In this paper, we examine the SFHs of NGC 1846 
and NGC 1783 using the same CMD reconstruction 
method previously applied to NGC 419 and NGC 1751 
(|Rubele, Kerber fc Girardll20ld ; iRubele et al.ll201lD . This 
method is significantly different from other analyses 



of the same clusters that ap p eared in the lit e rature 
dMackev fc Brobv Nielsenl 120071; IMucciarelli et al.l 120081; 



the r ecent observations of the open cluster Tr 20 by iPlatais et alj 
l2012|y 

2 The tidal radius and concentration index of NGC 1783 are not 
well constrained from the ACS/WFC data alone, due to its large 
radius. Its large r c , instead, is well constrained and evident even 
from a simple visual ins p ection of the HST ACS/WFC images 
used bv lGoudfrooii et all J2009I) . We note that there are two pre- 
vious fits of King profile to NGC 1783 in the literature, both pro- 
viding smaller values of r c : lElsonl l)l992T) finds r c = 4.9 pc = 20" 
using ground-based data, while IMucciarelli, Origlia fc Ferrarol 
[|2007l ) find r c = 24.5" from the shallower ACS/WFC images 
from SNAP 9891 (PI: G. Gi lmore). These small values p r obabl y 
explain why, in figure 3 of IKeller. Mackev fc Da Costal fcoill) , 
NGC 1846 and NGC 1783 appear as if they had the same r c , 
which is not correct. 



Mackev et all |200S| ; iMilone et al.M2009l : IGoudfrooii et ail 
20091 , l2011ak which concentrate on the MSTO region of 
the CMD. The CMD reconstruction method finds the best- 
fitting model for the entire CMD above a given magnitude 
cut, without giving particular weight to any subset of the 
observed stars. As such, it may be more affected by errors 
in evaluating the contribution from field stars or the pho- 
tometric errors of faint stars - aspects that, in any case, 
are taken into account in the method. On the other hand, 
by using the entire available data set, the CMD reconstruc- 
tion method potentially enhances the statistical significance 
of the detected SFH features. This is important considering 
that the number of stars along the MMSTOs is anyway lim- 
ited to a few hundreds, even for the most populous LMC 
clusters. Moreover, this aspect might be crucial in the anal- 
ysis of less populous clusters. 

In the following, we will use the excellent imaging and 
photometry of NGC 1846 and NGC 1783 available from 
HST/ ACS, together with new sets of evolutionary tracks 
and isochrones (Sect. [2J. Sect [3] will apply the CMD re- 
construction method of SFH-recovery to the surrounding 
NGC 1846 and NGC 1783 LMC fields (Sect.EH and cluster 
areas (Sect. [3T3|) . Sect. [4] discusses the results in the frame- 
work of present scenarios for the formation of MMSTOs. 



2 THE DATA 

2.1 Cluster imaging and photometry 

The data set used in this paper comes from GO-10595 
(PI: Goudfrooij), and consists of one short and two 
long exposures in F435W, F555W, and F814W with 
small dither patterns to avoid the gap between two 
ACS/WFC chips. A detailed des cription of the o b serva - 
tion s and photometry is give n in IGoudfrooii et all (|2009h 
and IGoudfrooii et al. (|2011bl ). Nevertheless, in this paper 
we use th e simultaneous ePSF fi tting technique as it de- 
scribed in lAnderson et al.l (|2008l 'l , which fits the PSF si- 
multaneously on all exposures / observations of the clus- 
ter. Differently from IGoudfrooii et all l|2009h . the Charge 
Transfer Efficiency ( CTE) correction was performed using 
iRiess fc Mackl |2004l ) formula (ACS-ISR 2005). The derived 
photomet ry was calibrated i nto th e Vegamag system as de- 
scribed in IGoudfrooii et all (|2009l ). 

The left panels in Figs. [1] and [2] show the spatial repre- 
sentation of the stars we analyse in this work, for the two 
clusters. The right panels show how the stellar density varies 
as a function of radius from the NGC 1846 and NGC 1783 
centres, taking into account only the stars of F814W < 22, 
for which the photometry is nearly complete. Based on these 
figures, we define Centre and Ring regions having radii and 
areas as tabulated in Table [T] For both clusters, the Centre 
regions have radii very close to the measured core radii r c . 
The Ring external radii are selected to be three times the 
Center radii (or about 3 x r c ), which include at least twice 
the number of stars as in the Centre regions. The figures 
also indicate the flattening of the stellar density for radii 
r > 2700 pix, which probably represents the regions which 
start being dominated by LMC field rather than by clus- 
ter stars. These radii were chosen as the inner boundary of 
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Figure 1. Left panel: Map of the stars used in this work, in the xy plane of the ACS/WFC images of NGC 1846. The scale is of 
about 0.05"/pfx. The observed stars have been grouped in areas corresponding to the LMC field (red) and, for NGC 1846, an inner 
"Centre" (green) and outer "Ring" (blue). Right panel: The logarithm of stellar density as a function of radius from the NGC 1846 
centre. Error bars are the random errors. In the top right corner we indicate the total numbers of stars used to build the profile, selected 
at F814W < 22 mag. 
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Figure 2. The same as Fig.Q] but for NGC 1783. 



the LMC Field for each cluster. Although these Field radii 
are poorly defined - especially for NGC 1783 - their stel- 
lar densities are clearly very small compared to the Centre 
and Ring regions. The field stellar densities are also similar, 
as expected for clusters located in the same portion of the 
LMC. 

Fig. shows the ACS data for the different regions of 
NGC 1846 and NGC 1783, in the F814W vs. F435W-F814W 
and F814W vs. F555W-F814W CMDs. These plots will be 
used as a reference in our analysis. 

The CMDs for the clusters show very clearly the broad 
main sequence turn-off, the composite structure of the red 
clump, and other well-known CMD features such as the se- 



quence of binaries parallel to the main sequence, and the 
RGB, subgiants, and early-AGB bump. A simple compari- 
son between the CMDs for the Centre and Field reveals that 
the field contamination in the Centre of both clusters is close 
to negligible. Indeed, the stellar densities in the Fields are 
about 25 times smaller than in the cluster Centres (see last 
column of Table[T} . Notwithstanding, it appears evident that 
the Field region in NGC 1783 is marked by a population with 
about the same turn-off as the cluster, which is probably in- 
dicating that areas dominated by the "pure LMC field" have 
not been reached in this case. We will evaluate the impact 
of this possibility further down in our analysis. 
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Table 1. Selected regions. The last 2 columns refer to stars with F814W <22 mag. 



Region 


radii 




area 


# stars 


mean density 




(") 


(pc) 


(arcmin 2 ) 




(arcsec -2 ) 


NGC 1846 Centre 


r< 20 


r< 4.8 


0.380 


2081 


1.52 


NGC 1846 Ring 


20 < r < 60 


4.8 < r< 15.4 


2.806 


5321 


0.527 


NGC 1846 Field 


r> 133 


r> 32.2 


2.454 


531 


0.060 


NGC 1783 Centre 


r< 33 


r< 8.0 


1.141 


6685 


1.63 


NGC 1783 Ring 


33 < r < 107 


8.0 < r< 25.9 


6.742 


6901 


0.284 


NGC 1783 Field 


r> 133 


r> 32.2 


2.212 


498 


0.063 





Figure 3. CMDs for NGC 1846 and NGC 1783. The selected regions are the same illustrated in Figs. Q] and [2] 



2.2 Assessing photometric errors and 
completeness 

In order to characterize the errors in the photometry and 
the completeness of the sample, we have performed a series 
of artificial star tests (AST) on the reduced images (see e.g. 
iGallart et alii 19991 ; [Harris fc Zaritsky||200 ll ). The procedure 
consists of adding stars of known magnitude and colour at 
random places in each exposure, and redoing the photome- 
try exactly in the same way as described in Sect. 12.11 The 



artificial stars are considered to be recovered if the input 
and output positions are closer than 0.5 pixels, and flux 
differences are less than 0.5 mag. In order to avoid the in- 
troduction of additional crowding in the images, artificial 
stars are positioned at distances much higher than their 
PSF width. So, our AST are distributed on a grid spaced 
by 20 pix, which is each time randomly displaced over each 
set of exposures. Importantly, the AST tests are repeated 
many more times in the central cluster regions, in numbers 
which are proportional to the density of stars brighter than 
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F814W < 22.5 mag. In this way, we have a better descrip- 
tion of the errors in the most crowded cluster regions, and 
are able to accurately describe their decrease with the radial 
distance from the centre. 

A total of 5.2 million ASTs were performed, with colors 
and magnitudes covering in an almost uniform way the CMD 
area of the observed stars and of the "partial models" to 
be used in the SFH analysis (see Sect. [3711 below). For the 
cluster centres, the 90 % completeness limit turns out to 
located at F814W~24.5, which is well below the position of 
the MMSTOs in both NGC 1846 and NGC 1783, as can be 
seen in Fig. |3] 

2.3 Stellar models 

An additional goal of this paper is to employ a new set 
of evolutionary tracks derived from the PAdova & tRi- 
este Stella r Evolution Code (PA RSEC) and extensively de- 
scribed in iBressan et all (|2012l ). They include updated in- 
put physics (opacities, equation of state, neutrino losses, 
etc.), and revised prescriptions for the convective processes, 
including microscopic diffusion in low-mass stars and an 
excellent description of helioseismic data. These improve- 
ments are expected to provide a more detailed description 
of CMD features and a more robust age scale than previ- 
ous versions of Pad ova tracks. In the mass interval of in- 
terest for this work, overshooting is assumed to operate 
with an efficiency of A c = 0.5 pressure scale heights (cf. 
IBressan. Chiosi fc Berteliilll98ll ). 

The init i al ch emical composition is derived from the 
ICaffau et all (|201ll ) new solar composition after enhancing 
the abundances of a elements by +0.2 dex. This mild a 
enhancement has just a minor effect in the shape of evo- 
lutionary tracks and isochrones, but it is necessary to re- 
produce the c hemical composition o f AGB stars observed in 
NGC 1846 bv lLebzelter et al.l (|2008l ). as will be discussed in 
Marigo et al. (in prep.). 

The stellar evolutionary tracks are transformed into 
isochrones and converted to the ACS/WFC Vegamag pho- 
tometric system u sing the transformations described in 
iGirardi et ail (|2008l ). 



3 RECOVERING THE SFH 
3.1 Overview of the method 

To recover the SFH from the ACS data, we use the 
same method applied to the SM C star cluster NGC 419 
ijRubele. Kerber fc Girard i 2010) and later impr oved with 
the LMC cluster N GC 1751 dRubele et all l201ll) : We use 
the StarFISH code (|Harris fc Zaritsky||200ll . |2004| ) to look 
for the linear combination of "stellar partial models" (SPM) 
that best-fits the Hess diagra of the observations, via the 
minimization of a x 2 -like statistics fcf. lDolphinl l2002l). The 
coefficients of this best-fit linear combination directly trans- 
lated into the SFH, and can be plotted in different ways - 
for instance, as the star formation rate as a function of age, 
SFR(t), plus the age-met allicity relation (AMR), [Fe/H](t). 



3 The Hess diagram is simply a representation of the stellar den- 
sity - number of stars per color-magnitude bin - across the CMD. 



SPMs are the basic building blocks in the method. 
They are theoretically-derived Hess diagrams of simple stel- 
lar populations spanning very small ranges of age and metal- 
licity. They are initially produced in a purely theoretical way, 
with the aid of the TRILEGAL population synthesis code 
(jGirardi et al.ll2005l ) and the PARSEC stellar evolutionary 
tracks. Then, these "perfect" SPMs are displaced by the dis- 
tance modulus and reddening to be tested, and degraded us- 
ing the distributions of incompleteness and photometric er- 
rors as derived from the ASTs. Specifically, all ASTs falling 
in the spatial region under consideration, and for each small 
box in the CMD, are grouped together and used to derive 
the two-dimensional error and completeness distributions, 
that then are used to blur the same boxes in the theoretical 
S PMs. Examples of this procedure are presented in figure 4 
of lRubele. Kerber fc Girardi (|2010l ) . 

In the case of star clusters, our customized version of 
StarFISH performs the following steps: 

(i) For each set of ACS frames, we first recover the best- 
fitting SFH of the Field region, exploring different values of 
Ay and (m — M)o, as described in Sect. 13.21 

(ii) From this best fitting solution and our set of stellar 
models, we generate a Field Stellar Partial Model (FSPM). 
This special partial model is (a) scaled to the cluster area to 
the analysed, and (b) degraded using the ASTs performed 
in the cluster area. Finally, this FSPM is included as a fixed 
component during the SFH-recovery of the Center and Ring 
regions. In this way, the SFH-recovery of the cluster area 
includes the best possible estimate for the LMC field con- 
tamination. 

(iii) We perform the SFH-recovery of the cluster area ex- 
ploring a wide range of extinction, distant modulus, and 
met allicity values (Sect. [TPll , 

(iv) This process is initially performed assuming a single 
value of binary fraction / (for binaries with mass ratios in 
the range between 0.7 a nd 1.0), and a fix ed present-day mass 
function (PDMF) from lChabrierl l|200ll ). Variations in these 
parameters are later explored. The default value of / is 0.3 
fo r the LMC field, an d 0.2 for the star clusters (cf. sect. 4.1 
in lRubele et al.ll201ll ). 

For the specific case of our NGC 1846 and NGC 1783 
data, we will use the entire CMD regions above F814W < 
22.5 mag. This ensures that we will be dealing with near- 
complete CMDs including all main sequence turn-offs up to 
the oldest possible ages. 

3.2 The SFH in the LMC Fields 

Deriving the SFH in the LMC fields is not a main goal in 
this paper, but it is both an interesting side-product of the 
present analysis, and a necessary step to reduce systematic 
errors in the cluster SFH. 

For both the NGC 1846 and NGC 1783 Fields, we de- 
rive the SFH usin g the same set of SPMs as defined in 
iRubele et all (|201ll ). We then run StarFISH to find the best- 
fitting solution to the observed CMDs, for a given value of 
distance modulus (m — M)o and extinction Av ■ Both F814W 
vs. F435W-F814W and F814W vs. F555W-F814W Hess di- 
agrams are used simultaneously in the process of \ 2 mini- 
mization. 

Figure [4] shows the map of Xmin ~ that is, the \ 2 value 
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Figure 4. \ 2 map for the Field best-fitting solutions, as a function of distance modulus and V-band extinction. The dashed and 
continuous black line show the 68 % and 95 % confidence levels for the overall best-fitting solution. The best-fitting model is located 
at (rn-M) = 1 8.685, Ay = 0.33 for the NGC f846 Field (left panel), and at {m-M) = 18.60, A v = 0.225 for the NGC 1783 Field 
(right panel). 




log(T/yr) log(T/yr) 

Figure 5. The SFH derived for the NGC 1846 and NGC 1783 Fields (left and right panels, respectively). In both cases, the top 
panel presents the best-fitting SFR(t) for the field (blue histogram), together with the random errors (lcr blue bar) and systematic errors 
(gray shadow) as estimated from the entire 68 % confidence level interval in the (m — M)o, Ay plane. The bottom panels present the 
mean age-metallicity relation (red and black point) with stochastic errors (red and black) and systematic errors (shaded regions). The 
green points show the center of distributions of the SPMs used to derive the SFH. 



to which StarFISH converges - for the solutions in the 
(m — M)o x Ay plane. We explored a range in both param- 
eters just extended enough to allow a clear identification of 
the absolute minimum and most of its 68 % confidence level 
interval. There are a couple of remarkable aspects in our 
results. 

First, the best-fitting models for the NGC 1846 and 
NGC 1783 Fields (left and right panels of Fig. g] respec- 
tively) present quite a different mean extinction, with Av = 
0.33 mag and Av = 0.225 mag, respectively. In comparison, 
for a radius of 12 arcmin a r ound the clusters NGC 1846 and 
NGC 1783, IZaritskv et al.l (|2004l ) derive extinction values of 
(A v ) = 0.59 ± 0.40 and (A v ) = 0.41 ± 0.42, respectively 
{{A v ) = 0.48 ± 0.32 and (Ay) = 0.31 ± 0.26 for cool stars) . 
The reddening maps by lHaschke. Grebel fc Duffaul (|201ll ) 
cover only the region of NGC 1846, providing (Ev-i) ^ 
0.06±0.075, which translates into (Ay) > 0.14 mag. Consid- 
ering the errors and the large dispersion in extinction values, 
all these values are consistent with each other. 

Second, the best-fitting distances for both fields just 
marginally agree one with each other, considering their 



68 % confidence level: indeed the distance modulus of 
(m— M)o ~ 18.625 mag represents, at the same time, a 
lower limit to the distance of the NGC 1846 field, and 
an upper limit to the distance of the NGC 1783 one. Al- 
though these distance measurements are perfectly compat- 
ible when we consider their 95 % confidence levels, we 
cannot refrain from noticing this unexpected result. Both 
fields have similar location in the Northwest portion of the 
LMC, relatively close to the line of nodes of the LMC disk. 
So, acco rding to recent results for the LMC disk geome- 
try (e.g. van der Marel fc Cionil 2001 1_ van der Marel et al.l 

m r 



l200a iNikolaev et alj|2004 iRubele et alj|2012h 7 they would 
be expected to have the same distance of the LM C center 
(of (m—M)o ~ 18.46 mag, see iRipepi et a"Hl2012l and refer- 
ences therein). Although the relative distances between the 
cluster fields and the LMC centre could be affected by sys- 
tematic errors, the relative distances between both cluster 
fields should be quite solid. 

We note however that both the distances and extinc- 
tion values for the field populations are consistent, within 
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their la uncertainties, with the values found for the cluster 
centres, as detailed in the next subsection. 

Finally, Fig. [S] presents the SFR(f) and AMR corre- 
sponding to the best-fitting solutions for both Fields. It is 
remarkable that the SFR(t) recovered for the NGC 1846 
Field presents features that are consistent with those 



dOlsen 


19991; Holtzman et al. 


19991; Smcckcr-Hane et al.1 


12002); 


Javiel. Santiago & Kerberl 120051) and ground-based 


data 


Harris & Zaritskvl 200U 


20091; iRubele et all l2012h. 



There is an initial period of star formation at ~ 10 Gyr, 
followed by a minimum at log(i/yr) ~ 9.7, and then a more 
extended period of star formation for ages between about 1 
and 4 Gyr. For younger ages, our data includes too small 
an area to set stringent constraints on the SFH; however, 
there are hints for significant SFR at ages of about 300 Myr 
(log(£/yr) = 8.5), and a well detected burst of formation of 
stars with ~ 10 Myr (log(t/yr) = 7.0). Error bars are too 
large to allow a meaningful quantitative co mparison with 
the re sults from lHarris fc Zaritskvl |2009T ) and IRubele et all 
(|2012h . for nearby regions of the IMC. 

The SFR(t) for the NGC 1783 Field, instead, presents a 
marked "burst" at ages between 1 and 2 Gyr, which, as we 
will see later, coincides with the ages of cluster formation. 
This is a further evidence that in NGC 1783 Field we are 
still sampling the cluster population. 

Concerning the AMR, the results for the NGC 1846 
and NGC 1783 fields are quite similar and con- 
sistent with those derived from IMC star clusters 
(jKerber. Santiago fe Brocatoll2007l : lHarris fc Zaritskvll2009l ) 
and for the IMC field using differen t sets of data (e.g. 
ICarrera et ai1l2008l : IRubele et al1l2012t . 

3.3 The SFHs in NGC 1846 and NGC 1783 

For NGC 1846 and NGC 1783 we proceed deriving the best- 
fitting SFH in the same way as performed for NGC 1751 
jRubele et al.ll201ll l In these cases, we assume that clus- 
ter stars present the same mean [M/H] value for all ages, 
since so far there is no evidence for si gnificant spreads 



in metallicity in such star clusters ( e.g. Mucciarelli et all 
l200Sl : IRubele. Kerber fc Girardl 1201(1 : IRubele et al.ll201lh . 

We have explored seven [M/H] values: —0.57, —0.54, —0.52, 
—0.49, -0.47, —0.44, -0.42. For each one of these mean 
[M/H] values, a box-shaped metallicity distribution is as- 
sumed, with a total width of A [M/H] = 0.025 dex. This 
spread is similar to the separation between the mean 
[M/H] values, and contributes to produce results that vary 
smoothly as a function of metallicity. 

The age interval covered by our SPMs goes from 
log(f/yr) = 8.9 to 9.4, which is much wider than the inter- 
val suggested by the position of NGC 1846 and NGC 1783 
MMSTOs. We initially adopt an age resolution (bin width of 
SPMs) of Alogi = 0.05 dex. So, for each set of parameters, 
we have a total of 11 partial models - 10 for the cluster, plus 
the FSPM described in Sect. 13. ll - completely encompassing 
the age interval of interest. 



3.3.1 The SFH for the NGC 1846 Centre 

Complete maps of Xmin f° r the NGC 1846 Centre, as a func- 
tion of (m — M)o, Ay and metallicity, are presented in the 



first seven panels of Fig. [6] These results are obtained for the 
age resolution of Alogi = 0.05 dex. We can notice that the 
best solutions are found in the metallicity interval between 
[M/H] = -0.52 and —0.42. The last panel shows the same 
kind of map for the Ring, but limited to the metallicity that 
provides the best fit for the Centre, namely [M/H] = —0.49. 
This value is in excellent agreement with the one derived 
from t he Ca n triplet of cluster members bv lGrocholski et al.l 
(j2006h . 

The best solution for the Centre is for (m — M)o = 
18.57, Ay = 0.26, with a Xmin = 0.55. Such a small Xmin 
is already an indication of an excellent fit to the observa- 
tional data. This best-fitting solution and map of residuals 
are also presented in the Hess diagrams of Fig. [7] Finally, 
the best-fitting SFR(i) is shown in the upper left panel of 
Fig.0 

To evaluate the errors for all involved parameters, we 
first find the correspondence between the Xmin value for each 
model and its confidence level. This correspondence was es- 
timated simulating 100 synthetic CMDs generated with a 
number of stars equal to the observed CMD, using the best- 
fitting SFR(i) and its parameters as the input for the simu- 
lations. So, after recovering the SFH for this sample of syn- 
thetic CMDs, it was possible to build the Xmin distribution 
and to establish the relation between their values and the 
confidence level. 

In the Xmin maps of Fig. |6j we superimposed the 68 % 
and 95 % significance levels for all the solutions for the Cen- 
tre. Only for [M/H] values between —0.52 and —0.42 we 
have solutions within the 68 % significance level of the best 
solution. Based on this figure, we determine (m—M)o = 
18.57 ± 0.07 and A v = 0.26 ± 0.05 for the cluster Centre 
(with random errors at the 68 % significance level). 

The left panel of Fig.[8]shows the SFR(i) for the cluster 
Centre together with error bars, as derived for the initial 
age resolution of Alogi = 0.05 dex. The most basic feature 
in this plot is that the SFR(i) is clearly non-null for two 
age bins, spanning the log(i/yr) interval from 9.15 to 9.25 
(ages from 1.41 to 1.78 Gyr). This result is not only valid 
for the best fitting model, but also across the entire 68 % 
significance level volume of the Ay vs. (m — M)o and [M/H] 
diagrams. Solutions within this volume are used to define 
the range of systematic errors, which is also depicted in the 
figure. Finally, we note that the two bins of non-null star 
formation are found even if we adopt less restrictive limits 
for the random errors, i.e. if we plot all solutions inside the 
95 % significance level. 

However, the very small random errors for the A log t = 
0.05 dex solution clearly suggest that the data presents the 
potential for a more detailed determination of the age distri- 
bution. The middle panel of Fig. |8]shows the solution when 
we adopt 20 SPMs separated by A log(age) = 0.025 dex, as 
obtained with the same data and methods, but for a fixed 
value of (m-M)o = 18.57, A v = 0.26, [M/H] = -0.49. Now 
the solution is clearly non-null for 3 age bins, spanning the 
log(i/yr) interval from 9.175 to 9.25 (ages from 1.50 to 1.78 
Gyr). As a result of the smaller number of stars per bin, ran- 
dom errors are larger than in the previous A log t = 0.05 dex 
cas^S 



Note that systematic errors are computed only in the case of 
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Figure 6. The first 7 panels show the maps of the Xmin obtained from the SFH-recovery of the NGC 1846 Centre, as a function of 
(m — M)o and Ay, for several [M/H] values (from —0.57 to —0.42). The black lines delimit the regions within the 68 % (dashed lines) 
and 95 % confidence levels (continuous lines) of the absolute best solution, which is found at [M/H] = —0.49 and is shown in the fourth 
panel. The last panel shows the same for the NGC 1846 Ring, but only for the metallicity of [M/H] = —0.49. The minimum Xmin are °f 
0.55 and 1.1, respectively, for the Centre and Ring. In both cases, the best-fitting solutions are found at (m — M)o = 18.57, Ay = 0.26. 



We proceed with this experiment, determining the 
SFH for an even better age resolution, namely A log t = 
0.015 dex. The results are in the right panel of Fig. [S] This 
time, the random errors are really large, as a consequence 
of the quite small number of stars per age bin. Under these 
conditions, the result of the SFH analysis could depend very 
much on the particular choice of limits for the age bins. In 
order to appreciate the SFR(f) for several bin positions (and 
always with the same resolution, of A \ogt = 0.015 dex), the 



the default age resolution of A log t = 0.05 dex, for which wc fully 
explored the possible interval of (m — M)q, Ay , [M/H], The basic 
reasons for not recomputing the systematic errors when adopting 
a better age resolution, are essentially: (1) The large CPU times 
needed to explore the entire (m — M)o, Ay , [M/H] interval. (2) 
Test runs of our software indicate that the best-fitting values of 
these parameters depend little on the age resolution A log t be- 
ing adopted. Indeed, they are strongly constrained by the CMD 
portions corresponding to the lower main sequence, RGB, and 
main body of the red clump, which arc equally well fitted at 
all age resolutions. Therefore, the systematic errors obtained at 
A log t = 0.05 dex shall be considered as indicative of those ex- 
pected at all age resolutions. 



age bins have been shifted progressively by steps equal to 1/5 
of the total bin width, that is, by 0.003 dex. The right panel 
of Fig. [S] is the result of plotting all these results together as 
a function of the bin central age. The SFR(f) appears much 
more continuous than in previous cases, and presents clear 
indications about two points: 

(i) The SFR(t) is non-null in the complete age interval 
from log(i/yr) = 9.18 to 9.25 (ages from 1.51 to 1.78 Gyr), 
which is in perfect agreement with the interval revealed by 
the SFR(£) of intermediate resolution, Alogf = 0.025 dex. 

(ii) There is a marked minimum in the SFR(t) for the 
age bins with log(t/yr) going from 9.22 to 9.23 (ages 1.66 
to 1.70 Gyr). Despite the large error bars, this minimum is 
statistically significant. It indicates, surprisingly, that there 
might have been a hiatus in the SFR(t) of NGC 1846, start- 
ing ~ 150 Myr after the first episode of star formation, and 
lasting for about 50 Myr. 

3.3.2 The SFH for the NGC 1846 Ring 

In the case of the NGC 1846 Ring we assume it has the same 
[M/H] as the cluster Centre, and explore the solutions in the 
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Figure 7. The Hess diagrams for the NGC 1846 Centre region (top panels) and Ring region (bottom panels). From left to right, the 
panels show the cluster F435W— F814W vs. F814W diagram, its best-fitting solution model, and the chi2 map. The same is done for the 
F555W-F814W vs. F814W diagram. 



(m— M)o vs. Ay plane, as shown in the last panel of Fig. [6] 
The assumption of the same [M/H] is natural for stars be- 
longing to the same (presumably chemically-homogeneous) 
cluster. In this case, the cluster Centre is taken as the ref- 
erence because it suffers less from field contamination and 
differential reddening, which are two potential sources of sys- 
tematic errors in the determination of [M/H]. Moreover, the 
cluster Centre presents low Xmin values, which means very 
good overall solutions. 

The best solution for the Ring turns out to be located 
at the same (m — M)o as the cluster Centre, and at virtually 
the same Av, as shown in Fig. [5] Experiments of SFH recov- 
ery at varying age resolution, as performed for the Centre, 
were repeated in the Ring. Also these results turned out to 
be remarkably similar to the Centre ones, as revealed by the 
bottom row of Fig. [HI In particular, the hiatus in star forma- 
tion at log(£/yr) = 9.22 is clearly present also in the Ring. 
We will further comment on this in Sect. [3] 



3.3.3 The binary fraction and mass function in 
NGC 1846 

Apart from the already-mentioned assumptions (constant 
[M/H], no differential reddening, etc.) our analysis also 
adopted SPMs built for a fixed value for the PDMF slope 
and fraction of unresolved binaries. The question arises 
whether our results can help to better constrain these pa- 
rameters, and their possible radial variations within the clus- 
ters. We test this by varying these parameters while keeping 
the same (m— M)o, Ay, and [M/H] as in the best-fitting 




log(M) [uj 



Figure 9. PDMFs used to derive the best-fit solutions for the 
clusters. The red, blue and green lines show the PDMF with 8 
different s lopes, a, from 0.4 to 3.6. The dot-dashed line shows the 
IChabriej 1120011 ) PDMF previously used to recover the detailed 
SFHs for the clusters and their fields. The dashed vertical lines 
indicate the mass interval relevant to this work. 



solution. In order to be more sensible to differences in the 
PDMF slope, we now cut the CMDs at a limiting magnitude 
of F814W< 24.25 mag, which corresponds to main sequence 
stars with ~ 0.7 Mq. Needless to say, by going to deeper 
magnitudes we are also including CMD regions of smaller 
completeness and with higher photometric errors. Although 
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Figure 8. Top row (a): The SFR(t) normalized to the average SFR over the complete age interval, for the NGC 1846 Centre, and 
derived for three different widths of the age bins, namely A log t = 0.05, 0.025 and 0.015 (from left to right, respectively). In the leftmost 
panel, the red histogram shows the SFR(t) for the best-fitting solution, while the orange bars show the random errors at the 68 % 
confidence level. The black dashed lines show the systematic errors inside the area of 68 % confidence level shown in Fig. [6] The central 
panel shows the same for a resolution twice as better, resulting on a more detailed SFR(i), but with random errors about twice as large. 
Finally, the rightmost panel shows the SFR(t) computed with a fixed resolution of A logt = 0.015, but moving the age bin centers at 
small steps of 0.003 dex in logt, so that the SFR(f) appears more continuous — and indeed, we change the SFR(f) plot to use continuous 
lines instead of histograms, in order to better illustrate these particular results. The continuous red line is the SFR(t) while the dashed 
orange lines illustrate the interval of random errors, at the 68 % confidence level. Bottom row (b): The same for the NGC 1846 Ring 
region. To provide a better comparison between the error bars among the different cases, the smaller sub-panels at the bottom of each 
panel show the errors (random, and also systematic in the case of A logt = 0.05 dex) as referred to the mean SFR(f) line. 



these processes are properly modelled by our ASTs, it is 
also true that the PDMF determinations we are doing here 
shall be considered of a more exploratory nature, than our 
previous determinations of the cluster SFHs. 

The binary fraction is parametrized by the fraction / 
of detached binaries with a mass ratio in the interval from 
0.7 to 1. These binaries cause the well-known sequence par- 
allel to the main sequence in the CMD, while binaries with 
a smaller mass ratio leave hardly any signature i n optical 
CMD s. The PDMF slope is parametrized by the ISalpeteJ 
(1955) slope, a (with dN/dM oc M a ), which changes the 
number ratio between the stars at and above the MSTO, 
and the fainter main sequence. PDMF slopes tested in this 
work are illustrated in Fig. [9] As a reference, the main se- 
quence stars at F814W ~ 24.25 mag, have ~07M a . For the 
mass range of interest here, the lChabrierl (|200lh PDMF has 
a slope of q — 2.2. 

In this exercise, we start from the best fitting solutions 
for an age resolution of 0.05 dex, and run StarFISH for all 



/ and a combinations, for a fixed distance and reddening. 
The results are illustrated in the Xmin map of Fig. II Ot sep- 
arately for the NGC 1846 Centre and Ring. In this figure, 
the Xmin maps were normalized to the minimum value of 
Xmin, to allow a better comparison between our results. The 
best-fitting solutions are found for about the same values 
of a and / for both Centre and Ring: (a = 1.2, / = 0.15) 
and (a = 1.6, / = 0.13), respectively. Since we do not per- 
form the error analysis, it is impossible to tell whether these 
small differences between Centre and Ring are statistically 
significant or not. Anyway, these values seem to indicate 
flatter PDMFs than the ones commonly found in galac- 
tic or extragalactic stellar c lusters, w h ich ty pically present 
PDMF s l opes close to the ISalpeteij (|l955l . 



Chabrierl |200ll, a ~ 2.2) one (see also e.g. |K 



20021 ; iBastian. Covey fc MeverjbOloT l . 



2. 35) or 
roupal |2001| . 



© 2012 RAS, MNRAS OOO.ITIUSI 



The SFH of NGC 1846 and NGC 1783 11 




Figure 10. Xmi„ map normalized to the minimum Xmin ^ or * ne NGC 1846 Centre (left panel) and the Ring (right panel) solutions, as 
a function of binary fraction / and PDMF slope a. 
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Figure 11. The same as Fig. IBI but for NGC 1783. The minimum Xmin ' s °f 1>39 and 1.45 for Centre and Ring, respectively. 



3.3.4 The SFH of NGC 1783 Centre and Ring 

We proceed in the analysis of NGC 1783 exactly in the same 
way as for NGC 1846. Complete maps of Xmin f° r the Cen- 
tre, as a function of (m— M)o, Ay and metallicity, are pre- 
sented in the first 7 panels of Fig. 1111 The best solution 
is found again for [M/H] = -0.49, with (ra-M)o = 18.57, 
Av = 0.22, and a Xmin = 1-39. Solutions within the 68% con- 
fidence level span the range of metallicities between —0.47 
and -0.54. 



The metallicity of [M/H] = —0.49 turns out to coincide 
with the one previously determined for NGC 1846, which 
is not unexpect ed given their same ag es. High resolution 
spectroscopy bv lMucciarelli et al.l (|200&t ) instead indicates a 
value of [Fe/H] = -0.35 ± 0.02. 

Based on Fig.Ql] we determine (m-M) = 18.57±0.07 
and Ay = 0.22 ± 0.05 for the cluster Centre (with random 
errors at the 68 % significance level). 

Using [M/H] = —0.49 for the Ring, we find good so- 
lutions in quite a similar region of the (m — M)o vs. Ay 
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Figure 12. The same as Fig. but for NGC 1783. 
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plane as for the Centre (see last panel in Fig. The 
overall best solution is at (m — M)o = 18.56, Ay = 0.20, 



and has a 



1.45. We however take the solution at 



(m — M)o = 18.57 as the reference one, which presents a 
minimum Xmm a t Ay = 0.22 just as for the Centre. The 
best-fitting solutions and map of residuals are also presented 
in the Hess diagrams of Fig. 1121 while the SFR(t) are illus- 
trated in the left panels of Fig. 1131 

The panels of Fig. [13] show the SFR(t) for the cluster 
Centre and Ring together with error bars, again at several 
different resolutions A log t. These solutions are remarkably 
similar to those already found for NGC 1846, in several as- 
pects. Age intervals of non-null SFR(t) are essentially the 
same. In particular, also in NGC 1783 we find indication 
for a hiatus in the SFR(t) taking place ~ 150 Myr after 
the initial burst of star formation. Comparison between the 
middle panels of Fig. 1131 however, reveals something new: 
the SFR(t) tends to be slightly older in the cluster Centre 
than in the Ring. In particular, the Ring seems to present a 
less marked peak of star formation at older ages. 

To verify whether this difference could be due to an in- 
correct account of the field contamination in the NGC 1783 
Ring region, we make an additional test, illustrated in 
Fig. 1141 The contribution of the field is artificially in- 
creased/decreased by 50 %, with respect to the contribution 
expected for the Ring area, and then the SFH recovery is re- 
peated. The results are almost indistinguishable from those 
in the bottom-right panel of Fig. 1131 indicating that this is 
not the factor driving the different results between Center 
and Ring in NGC 1783. 



3.3.5 The binary fraction and mass function in 
NGC 1783 

As for NGC 1846, also for NGC 1783 we take the best-fitting 
solutions for Centre and Ring and vary its binary fraction 
/ and PDMF slope a. The maps of Xmin are presented in 
Fig. 1151 The results this time appear noticeably different 
between Centre and Ring: While for the Centre the best- 
fitting values are a — 1.3 and / = 0.20, for the Ring they 
are a = 2.0 and / = 0.15. These values of a provide hints 
of a possible mass segregation in this cluster. 

The interpretation of this finding is not straightforward. 
Due to CPU time limitations we did not perform a thorough 
exploration of the parameter space, neither the detailed 
analysis of random and systematic errors that would be re- 
quired to define the significance of the detected PDMF vari- 
ation. Anyway, we think the results are interesting enough 
to be mentioned here, considering the novelty of determin- 
ing a in different regions of clusters with MMSTOs. Since 
we find a smaller contribution of the older stars in the Ring, 
the deficit of less massive stars in the Ring is more likely 
associated to the younger population. 



4 DISCUSSION 

Our analysis produces a number of results, that we now com- 
pare to previous results for the same clusters, and frame in 
the present scenarios for the production of multiple popula- 
tions in Magellanic Cloud clusters. 
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Figure 13. The same as Fig.|8]but for NGC 1783. 




Figure 14. The high-resolution SFH derived for the NGC 1783 Ring, now assuming a field contamination decreased/increased by 50 % 
(left /right panels, respectively), with respect to the values expected for the Ring area. 



4.1 Comparing the SFH between NGC 1846 and 
NGC 1783 Centres 

As easily noticeable by comparing the top right panels of 
Figs. [S] and [H the timescales for the SFH in NGC 1846 and 
NGC 1783 Centres are surprisingly similar. They both seem 
to have formed stars for a total period of about 0.3 Gyr. 
They both started the star formation activity in a marked 
burst peaked at log(t/yr) = 9.24 or 9.25 (t = 1.78 Gyr). The 
FWHM of this marked peak is, in both cases, of 0.02 dex 



in \ogt, which is indistinguishable from the actual resolu- 
tion of Alogt = 0.015 dex of the SFH-recovery method. In 
both clusters, this strong peak is followed by a short hiatus 
in SFR(t), which lasts long enough to be detected by the 
method, that is, at least 0.02 dex in logt (or 70 Myr). Then 
the SFR(t) proceeds up to more recent times, becoming null 
at ages of log(t/yr) = 9.17 (1.48 Gyr). In the following, for 
obvious reasons we will refer to the marked peak observed in 
the SFR(t) as the first generation, and the extended period 
of star formation, after the hiatus, as the second generation. 
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Figure 15. The same as Fig.[[0] but for NGC 1783. 



It should be clear however that more than two generations 
may be present, especially in the later periods of star for- 
mation (younger stellar ages). 

The presence of separated periods of star formation in 
these clusters has been already advanced by other authors, 
using essentially the same set of HST images, but differ- 
ent da t a reductions and met hods. iMackev fc Brobv Nielsen! 
|2007lY IMackev et ail |200Sl ') and IMilone et al.l (|2009T l no- 
ticed that the MSTO region in NGC 1846 was clearly split 
or bimodal, hence suggesti ng the presence of t wo m a in pop - 
ulations separated in age. iGoudfrooii et all (|2009l . l2011bT ) 
instead derive a continuous distribution of stars across the 
MSTO region, hence concluding that the SFR(t) likely 
proceeded continuously over a period of ~ 30 Myr. For 
NGC 1783. lMucciarelli. Origlia fe Ferrard J2007I) find n o ev- 
idence of a broad turn-o ff, while IMackev et al. |2008l ) and 
IGoudfrooii et al.l (|2011bh c laimed a MMSTO w ith a smooth 
age distribution, whereas IMilone et al.l (|2009T l suggests a 
double MSTO. Our paper, instead, would agree more with 
the interpretation of a split MSTO in both clusters. Notice 
that the different results cannot be attributed to the dif- 
ferent data reduction only. The main different between our 
analysis and the previous ones is in the fact that we use 
the information in the entire CMD, and not only the one 
across the MSTO, in order to derive the age distribution. 
Our method also attempts to fit other age-sensitive features 
like the subgiant branches and red clumps - very well drawn 
in the CMD of both clusters, see Fig. [jj] 

As already mentioned in Sect.fi] the most evident differ- 
ence between NGC 1846 and NGC 1783 is in their different 
radii, while most of their other parameters are quite similar 
- including the total mass, mean metallicity, and level of 
contamination by the LMC disk. Both clusters are among 
the most massive of their ages in the LMC. It is not hard 
to imagine that they formed at about the same age, as a 
result of the same sort of dynamical process (or the same 
large-scale dynamical event) in the LMC. All the above- 
mentioned aspects point to a common process operating in 
a similar way in both clusters, from the onset of star forma- 
tion up to later stages. 

Among the scenarios advanced to explain prolonged star 
formation in LMC star clusters, two are particularly worth 
of being commented here: First, from dynamical arguments 
and the inferred escape velocity a t an age of 10 My r of sev - 
eral clusters exhibiting MMSTOs, IGoudfrooii et all (|2011aT ) 
suggest that second-generation stars start being formed from 



cluster material shed by first-generation stars featuring slow 
stellar winds. In our results, the second generation appears 
after a period of time to be located between 160 and 220 Myr 
after the first generation. This time lag corresponds to the 
main-sequence lifetimes of stars with masses between 4.0 
and 3.4 Mq, which are expected to quietly end their nu- 
clear lives as mass-l osing AGB stars with low -velocity winds 
(~15-20 kms -1 ; lHabing fc Olofssonlliool ). Thus, our re- 
sults for the duration of the SFH hia tus are certainly com- 
patible with IGoudfrooii et all (|2011al Vs results. 

A somewhat simila r scenario has been advanced by 
IConrov fc Spergell |201ll ). who suggest that the supernovae 
of type II and the prompt type la from the first-generation 
stars, first clean up the interstellar medium of the cluster; 
the interstellar medium is then reformed at the cluster cen- 
tre as intermediate-mass stars start shedding their envelopes 
at ages larger than 100 Myr. The accumulated gas then 
cools as the Lyman- Werner photon flux from the cluster 
stars drops after a few 100 Myr, which allows molecular hy- 
drogen and stars to form. This starts the second-generation 
of star formation, after which no third generation can fol- 
low because late SN la from the first generation clean up 
the cluster gas again. Although detailed in the description 
of the pr ocesses that can a ll ow/in terrupt star formation in 
clusters, IConrov fc Spergell l|201ll Vs scenario do not make 
very precise predictions about the timescales of the differ- 
ent processes involved. Anyway, it is clear that our findings 
regarding the SFH, with the observation of first and second- 
generation stars, are also in agreement with their descrip- 
tion. In addition, our observations indicate that in these 
~ 10 Mq clusters the final dissolution of the central inter- 
stellar medium by late SN la takes place after 0.3 Gyr. In the 
SMC cluster NGC 419, using the same method we find evi- 
dence fo£_a_totalinterval_of_^t at least twice as 
long (jRubele. Kerber fc Girardll20ld ). These numbers may 
provide important constr aints for a further refinement of the 
IConrov fc Spergel (|201lh scenario. 



4.2 Comparing Centre and Ring SFHs 

Another finding from our methods regards the similari- 
ties/differences between the SFHs at the cluster Centres and 
Rings. As indicated by the right panels in Figs. 151 and !13l we 
find hints of a different SFH between Centre and Ring for 
NGC 1783, in which the Ring has a much less pronounced 
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first-generation peak. For NGC 1846, instead, there is no 
suggestion of a similar effect. 

This result i s partic ularly surprising since results from 
iGoudfrooii et all (|2011al 'l. derived from the same data, seem 
to suggest exactly the opposite: after drawing boxes in the 
CMD for these and other clusters, they counts stars in the 
"upper and lower MSTO regions" , finding that both kinds of 
stars present a very similar radial distribution in NGC 1783, 
while in NGC 1846 the upper MSTO stars are clearly more 
centrally concentrated. 

As already commente d, the present analysis uses all 
stars in the CMD, whereas IGoudfrooii et all l|2011al) explic- 
itly avoids dealing with background/field subtraction issues 
by excluding areas in the CMD where the background was 
found to contribute more than 20% to the star counts. Their 
goal is to d eal with quite robust star counts as function of ra- 
dius. Thus. IGoudfrooii et al.l (|2011al ) do not use a large part 
of the MS, and the top-left extreme of the MSTO, which 
instead are included in the present analysis. In addition, 
IGoudfrooii et al.l (|2011al Vs method is model-independent. 

Our approach instead aims to maximize the number 
statistics by including all stars irrespective of their origin, 
and, additionally, to limit the impact of subjective choices 
regarding the stars included in the analysis. This attempt 
cannot be free of subjective choices, however. For instance, 
we had to decide, quite subjectively, the extent of Centre 
and Ring regions. Also, we use a particular set of isochrones, 
under the implicit assumption that they accurately describe 
all relevant phases of stellar evolution. Inaccuracies in these 
isochrones would unavoidably appear as systematic errors 
in all of our results. 

These contrasting results clearly make somewhat un- 
certain any attempt to use the radial trend to discuss the 
scenario for the formation of second-generation in these clus- 
ters. Formally speaking, the results for NGC 1846 - no dif- 
ference between Centre and Ring SFH - seem just to indi- 
cate that we were not able, with our definition of Centre 
and Ring, to detect the expected concentration of the sec- 
ond generation towards the centre of this cluster. This is 
not surprising. What is surprising is the trend derived for 
NGC 1783, which would suggest that, in less concentrated 
clusters like this, second-generation stars are spatially more 
spread than first-generation ones. No scenario for the forma- 
tion of multiple populations in clusters seem to favour such 
a trend. This point clearly deserves more accurate analysis. 

As for the determination of the PDMF and binary frac- 
tions, the basic results is that we find hints of mass segrega- 
tion in NGC 1783, at a level which is certainly larger than 
in NGC 1846, but which, however, may not be statistically 
significant. NGC 1783 is also the cluster with the largest ra- 
dius, the larger concentration index c, and also the one for 
which Centre and Ring seem to present a different SFR(t). 
It is tempting to suggest a correlation between all these dif- 
ferent aspects, but any conclusion on this is hampered by 
the uncertainties in our determin ations of a and /, in th e 
determination of c for NGC 1783 (IGoudfrooii et al.ll2011bT ). 
and on the fact that we have analysed just two clusters. It 
would be very interesting to extend the same kind of deter- 
mination to a much larger sample. 

At first sight, our clusters seem to follow the anti- 
correlation between the concentration index c and the slope 
of the global mass function, which is detected in Galactic 



globu lar clusters (GGCs ; |Pe Marchi. Paresce fc Pulonej 
12007? ) and interpreted as the result of residual-gas 
expulsion in initially mas s- segreg ated star clusters 
|Marks. Kroupa fc Baumgardtl l200Sf ) : indeed, we find 
the steepest mass function (larger a) in NGC 1783, the 
cluster with the larger concentration index. We can also 
confirm that, after taking into account the different defi- 
nitions of c and a (i.e. the use of logarithm for c, and the 
multiplication by —1 in the case of a) by differe nt authours, 
we fin d values comparable to those derived by iGlatt et alj 
(|201 ll ) for 6 SMC intermediate-age and old clusters; these 
values place them from ~ 2 to 3 dex bellow the a vs. c 
relation typical of GGCs. The interpretation of these trends 
is not easy, since all these clusters have been analysed using 
different methods and heterogeneous dat a , and , moreover, 
the mean |Pe Marchi. Paresce fc Pulonel ^OOTm 's relation 
is usually derived for entire clusters, and not for different 
cluster regions as in our case. 



4.3 Concluding remarks 

Together with our previous findings for th e SMC clusters 
NGC 419 (iRubele. Kerber fc Girardll2010D and NGC 1751 



(|Rubele et al.ll201ll ) using the same techniques, the present 
results for NGC 1846 and NGC 1783 confirm that multi- 
ple episodes of star formation provide a good quantitative 
description of the observed CMD features. At this point, 
research in the field should concentrate on carefully age- 
dating the different populations, and in identifying features 
that can favour/disfavour the possible scenarios for their for- 
mation. In this sense, the present work has been partially 
successful. Without imposing any a priori limitation or pre- 
selected shape to the SFH, we have found the presence of 
separated first and second-generation episodes of SFH for 
both clusters. This observation is compatible with simple 
schemes where second-generation SFR(i) can only start af- 
ter a period of a few 100 Myr, when the prompt supernovae 
and t he Lyman- Werner flux f rom first-generation stars are 
over (|Conrov fc Spergell l201lh . The stellar ejecta accumu- 
lated to form the second-generation should come partially 
from 3.4-4.0 Mq stars of the first generation, and be shed 
through slow AGB winds, in accordance with the correlation 
between the pres ence of MMSTOs and v e 3C at an age with 
10 Myr found by IGoudfrooii et all (|2011al ). 

However, we find an intriguing central concentration of 
first-generation stars in the cluster NGC 1783, which is op- 
posite to what predicted by the above-mentioned scenarios. 
This cluster has a large core radius and concentration index, 
which may be associated with the observed gradients in the 
stellar population. The bad definition of the cluster field 
does not seem responsible for the observed trends, because 
the field contribution is anyway very small. However, clarifi- 
cation of this issue may requires the analysis of larger areas 
around the cluster, and of more clusters of similar age and 
different structural parameters. These steps will be pursued 
in subsequent papers. 
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